Simple proof of the robustness of Gaussian entanglement in bosonic noisy channels 
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The extremality of Gaussian states is exploited to show that Gaussian states are the most robust, among all 
possible bipartite continuous-variable states at fixed energy, against disentanglement due to noisy evolutions 
in Markovian Gaussian channels involving dissipation and thermal hopping. This proves a conjecture raised 
recently in [M. Allegra, R Giorda, and M. G. A. Paris, Phys. Rev. Lett. 105, 100503 (2010)], providing 
a rigorous validation of the conclusions of that work. The problem of identifying continuous variable states 
with maximum resilience to entanglement damping in more general bosonic open system dynamical evolutions, 
possibly including correlated noise and non-Markovian effects, remains open. 
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Continuous-variable (CV) systems such as light modes 
and ultracold atomic ensembles 1 1 1 provide advantageous re- 
sources to achieve unconditional implementations of quantum 
information processing 121, ranging from teleportation proto- 
cols 1 3 1 to quantum key distribution |4 1 and one-way quantum 
computation |5|. Gaussian states and Gaussian operations, 
that represent respectively the most natural and easily control- 
lable light states as well as the set of manipulations efficiently 
realizable by linear optics, have traditionally occupied a privi- 
leged role in all such implementations. Furthermore, by virtue 
of their mathematical simplicity compared to general states 
living in infinite-dimensional Hilbert spaces, bosonic Gaus- 
sian states have been and are the preferred testgrounds for a 
broad variety of investigations into the structure, nature, and 
dynamics of C V entanglement and quantum correlations 1 6 1 . 

Some prominent limitations of a Gaussian-only toolbox 
have been however recently exposed in several contexts |7|, 
stimulating vigorous theoretical and experimental research 
into the realm of non-Gaussian state engineering and charac- 
terization L8J, to assess and harness the potentially enhanced 
performance of de-Gaussified CV resources for quantum tele- 
portation f9l, entanglement distillation fTOl, parameter esti- 
mation 1 11 1, universal quantum computation El, non- 
locality tests 113, etc. Still, one of the most powerful fea- 
tures of Gaussian quantum states is their extremality 1 14] in 
the space of all CV states, that allows to formulate valuable 
bounds on suitable entanglement measures and entropic de- 
grees for a general non-Gaussian state g, based on the corre- 
sponding (easier to compute) properties of the Gaussian state 
a with the same first and second statistical moments as g. This 
has important consequences for the security of CV quantum 
key distribution ifTSl . 

The transmission of one-mode and multi-mode, possibly 
entangled beams between distant locations, a basic neces- 
sity for the realization of a distributed quantum communica- 
tion network |16|, is unavoidably affected by various types 
of noise. While phase diffusion (dephasing) dampens the co- 
herences in the Fock basis transforming Gaussian states into 
non-Gaussian ones, the exposition to dissipative losses and 
thermal hopping yields an instance of a Gaussian channel, that 
preserves the Gaussian character of the states sent through it. 

Recently, Allegra, Giorda and Paris have attempted a com- 



parison (by means of several separability criteria) of the deco- 
herence effects due to the latter types of channels, on special 
families of Gaussian as well as non-Gaussian two-mode en- 
tangled states, all initially having the same mean energy ifTTl . 
Among the main findings of their work, they formulate a con- 
jecture on the "robustness of Gaussian entanglement", accord- 
ing to which initial Gaussian pure states should be the ones 
whose entanglement is the last to vanish, compared to other 
classes of initial non-Gaussian pure states with equal energy, 
when they are all evolving through a noisy Markovian channel 
with loss and thermal hopping. 

In this Brief Report, we prove the conjecture of Ref. i fTTl 
by resorting to the extremality of Gaussian states 1 14], in par- 
ticular to the fact that Gaussian states maximize entanglement 
among all CV states with fixed energy 1 18]. Our simple re- 
sult establishes in full generality that no advantage can be 
gained by encoding information into entangled non-Gaussian 
states when they have to be transmitted through Gaussian 
noisy channels (under dissipation and/or amplification), and 
that Gaussian states are, hence, the most robust CV states in- 
deed against photon losses and thermal hopping. An outlook 
on the problem of identification of the most robust CV entan- 
gled states for more general bosonic noisy evolutions, includ- 
ing dephasing channels, is provided in detail. 

We consider general pure two-mode CV states ^12(0) with 
fixed initial mean energy no = Tr[^i2(0)(ajai + a2a2)], 
where aj , a j , are the ladder operators of mode j = 1,2 sat- 
isfying the canonical commutation relations [ai,a^] = dij. 
Among all the possible input states, the two-m ode squ eezed 
state ^^2(0) = Wii^l where = A^Vl - A2|n,n) 
and no = /(I — A^), is the only Gaussian instance (up to 
local unitary operations). We let each possible initial state un- 
dergo a dissipative evolution through the channel described by 
the following master equation 1 17], 

Qi2{t)= ^ ^NjL[a]]g,2{t)^^{N,^l)L[aj]gi2{t), 

(1) 

that encompasses losses and thermal hopping in the presence 
of nonclassically fluctuating local environments. The dot de- 
notes time-derivative and the Lindblad superoperator is de- 
fined as L[0]g = 20 gO^ - O^Og - gO^O. Here F is a 
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FIG. 1: (Color online) Scheme of the channel in Eq. ([T]). 



loss coefficient and Nj is the mean photon number in the sta- 
tionary (thermal) state of each mode. The master equation ([T]) 
admits the operator solution |[T7l|T3: ^12 (^) = A^^i2(0) = 
Tr34[[/t(^(0) ®i'M)Ut], where At denotes the evolution map 
corresponding to the noisy channel; 3, 4 are two additional fic- 
titious modes in a thermal state ^^34 = vs^v^'^Ut = Uis{(t)'S) 
^24(0) and Uij{Ct) = ex.p{Ctalaj - Qo}-ai) amounts to a 
beam-splitter transformation with Ct = arctan(e^^ — 1)^/^. 

The action of the noisy evolution in Eq. ([T]) can be thus sim- 
ply modeled as letting each mode interfere with an ancillary 
thermal mode via a beam-splitter with transmissivity cosC^, 
and eventually tracing over the ancillae, as depicted schemat- 
ically in Fig.[T] In this notation, A^i (A^2) is the mean photon 
number of the ancillary mode v^, (^4)- In the Heisenberg pic- 
ture, the beam-splitter transformation acts on the mode oper- 
ators as follows ||2Q1| 



; COS Ct 



sinO 



(2) 



where as, 4 are the annihilation operators of the ancillary 
modes. We can calculate the mean energy of the generic state 
^12 (t) at any time t during the evolution: 



Tr 



= no cos^ Ct + {Ni + N2) sin^ Ct , 



(3) 



where all the additional cross-terms vanish because the an- 
cillary thermal states have zero first moments, entailing 
(<^3,4)i/34 = (<^3 4)1/34 = 0- What is important here, is that 
the energy of the evolving states ^12 (^) at any time only de- 
pends on the initial energy and on the channel parameters, i.e., 
if a set of states ^12(0) enter the channel all with the same en- 
ergy, as it is the case in 117 1 and in the present discussion, their 
energies will remain equal to each other (yet globally affected 
by the open system dynamics) throughout the whole evolu- 
tion, regardless of the specific form of the states and of their 
Gaussian or non-Gaussian nature. Let us also remark once 
more that the channel in Eq. ([T]) is a Gaussian channel, i.e. the 
state gf2{'^) evolved from the two-mode squeezed state will 
remain Gaussian during the whole evolution. The final impor- 
tant ingredient we need is the extremality of Gaussian states. 



If one adopts a bona fide measure of entanglement £ II2TII , 
which is continuous and strongly superadditive, such as the 
distillable entanglement, or the squashed entanglement 1221 . 
then the Gaussian state (t) with energy fit is the most en- 
tangled among all possible CV states ^12 (^) with the same 
energy |[l4||T8l, as a simple consequence of the maximum en- 
tropy property of Gaussian states |23|. 

Collecting all the above observations, we can conclude that 
^[Qi2{^)\ ^ ^[Qi2{t)] at any time t, in other words the Gaus- 
sian two-mode squeezed state starts in advantage, in terms of 
entanglement, over any other CV state with the same energy, 
and maintains its advantage throughout the whole noisy evo- 
lution. Given that (as soon as A^i , A^2 7^ 0) the channel in 
Eq. ^ takes a finite "separation time" r[^i2(0)] to erase all 
the entanglement in any initial quantum state ^12(0) jlTl . the 
entanglement in the Gaussian state must be the last one to van- 
ish, i.e., r[^f2(0)] > ^[^12(0)] for all CV states ^12(0) with 
the same mean energy as (0) . This proves the robustness of 
Gaussian entanglement conjecture of Ref. 1 17 1, for continuous 
and strongly superadditive bipartite entanglement monotones, 
in the general case of Ni ^ N2 and without the need to re- 
strict ^12(0) to be in the form of a photon-number-entangled- 
state (i.e., a two-mode state whose one-mode marginal density 
matrices are diagonal in the Fock basis). Actually, being es- 
sentially a consequence of extremality, the present result holds 
for any set of, not necessarily pure, iso-energetic input states 
^(0) of an arbitrary number of modes that include a Gaussian 
one ^^(0), and for general Markovian Gaussian channels. 

All the speculations of Ref. ifTTl that were based on this 
conjecture can be now regarded as valid and truthful: Gaus- 
sian entanglement is extremal in terms of robustness against 
decoherence due to noise and dissipation. It is important to 
notice, however, that if entanglement is measured e.g. by the 
(logarithmic) negativity II2TII . which does not satisfy the hy- 
potheses of the theorem in fTT], then the extremality of Gaus- 
sian states can be violated by tiny amounts 1 14 , 18 1, hence nar- 
row regimes might be found where a particular non-Gaussian 
state could overcome ^12 (^) the entanglement hierarchy 
based on negativities. This is nonetheless to be regarded more 
as a pathological feature of the employed entanglement mea- 
sure, rather than as a physically meaningful violation or coun- 
terexample to the robustness conjecture. Let us remark once 
more that our general proof in fact holds for any continu- 
ous and strongly superadditive bipartite entanglement mono- 
tone (such as the squashed entanglement 1221 ). It becomes 
quite clear, henceforth, how the degree of non-Gaussianity - 
measured as distance from a Gaussian reference state with 
the same covariance matrix 1241 - of different instances of 
^i2(0)'s is set to play a key role in assessing their robustness 
to noise and their entanglement evolution, as evidenced by 
several examples in Ref. ITTl . 

Our analysis generalizes the one in Ref. fTSl, where the 
two-mode squeezed Gaussian state had already been identified 
as "the most robust state of light" against evolutions in purely 
dissipative channels [i.e. with A^^ = A^2 = in Eq. ([T])]. 
These results do not straightforwardly extend to more gen- 
eral scenarios in which other types of noisy channels are taken 
into account, that can violate Gaussianity and/or Markovian- 
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ity. In these cases, the Gaussian nature of transmitted states is 
not preserved, and the state evolved from ^^2(0) immediately 
loses its privileged role for t > 0, entering the competition 
among the ^i2(^)'s as a peer. This is the case, for instance, if 
purely local phase diffusive decoherence channels are consid- 
ered, which are non-Gaussian channels modeled by a master 
equation of the form 

^i2(t)= ^ ^L[a]aj]g,2{t). (4) 

For such noise models, the coherences in an initial two-mode 
squeezed state are degraded undergoing a non-Gaussian evo- 
lution, eventually driving the two modes into an uncorrelated 
product of thermal states, each corresponding to the reduced 
one-mode density matrix of the input two-mode state. In this 
scenario, there is no apparent hierarchy regulating the dynam- 
ical comparison between Gaussian and non-Gaussian initial 
states ^12(0), when energy or alternative resources (i.e., the 
squeezing degree) are fixed at time t = 0. More complicate is 
the situation when a two-mode dephasing channel, with cor- 
related noises affecting each mode, is taken into account; in 
that case, not even the marginal states are preserved. 

To formulate an outlook, we feel that even more challeng- 



ing, yet surely important, would be to investigate the case of 
more realistic channels where dissipation, thermal hopping, 
as well as phase diffusion and possibly memory effects and 
correlated noise may all simultaneously take place. Some of 
these mechanisms are currently investigated, for instance, in 
the context of unveiling quantum coherence effects and noise- 
assisted energy transport phenomena in biological systems 
125 1 . The problem of entanglement evolution and robustness 
turns thus into a highly non-trivial arena, and since the power- 
ful extremality results 1 14 1 are of little or no use at all in these 
more general cases, a numerical approach along the lines of 
Ref. 1 17 1 should be the preferable avenue to pursue. This will 
be the subject of further work. The main result of this Brief 
Report should then be recovered in the limit of negligible de- 
phasing, uncorrelated noises and memory less channels, but 
we feel that no reliable prediction can be made at this stage 
about which type of state, in the uncountable arena of CV 
states, might gain the crown of the most robust state of light 
in the intermediate regime where all the diverse noise effects 
are comparable. 
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